VERY SINGULAR SIMILARITY SOLUTIONS AND HERMITIAN 
SPECTRAL THEORY FOR SEMILINEAR ODD-ORDER PDES 



R.S. FERN ANDES AND V.A. GALAKTIONOV 

Abstract. Asymptotic large- and short-time behaviour of solutions of the linear dis- 
persion equation 

ut = Uxxx in K X R+, 
and its {2k + l)th-order extensions are studied. Such a refined scattering is based on a 
"Hermitian" spectral theory for a pair {B,B*} of non self-adjoint rescaled operators 

B = + ^ yDy + and the adjoint one B* = _ i yj^^ 
with the discrete spectrum (t(B) = ct(B*) = {A; = — -I, I = 0, 1,2, ... } and eigenfunc- 
tions {ipiiv) = '•^^ DyA\{y), I > 0}, where Ai{y) is Airy's classic function. Applications 
to very singular similarity solutions (VSSs) of the semilinear dispersion equation with 
absorption, 

us{x,t) ^t^^ f{j^) : Ut = Uxxx - lulP-'^u in M x M+, p > 1, 

and to its higher-order counterparts are presented. The goal is, by using various tech- 
niques, to show that there exists a countable sequence of critical exponents {pi = 
1 + ^ = 0, 1, 2, ...} such that, at each p = pi, a p-branch of VSSs bifurcates from the 
corresponding eigenfunction of the linear operator B above. 



1. Introduction: semilinear odd-order models, history, and results 

1.1. Basic dispersion models and applications. As a first basic model, we will study 
higher odd-order partial differential equations (PDEs) of the form 

(1.1) ut= {-l)''+^Df+^u + g{u) in M x M+, k = 1,2, ... , 

with bounded integrable initial data u{x,0) = uo{x) in M. Here D"^ = (d/dx)'^ denotes 
the m^^ partial derivative in x. The odd-order semilinear dispersion equation (11.11) can 
be considered as a counterpart of the better known semilinear higher-even-order parabolic 
PDE of reaction- diffusion type, 

(1.2) ut = {-lf~^^Dl''u + g{u) in M x M+, k = l,2,.... 
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For k = 1, (11.21) becomes a standard react ion- diffusion equation from combustion tlieory 



(1.3) Ut = Uxx + g{u) in M X M_|_, 

to wliicli dozens of well-known monographs are devoted to. These parabolic equations 
indeed belongs to an entirely different type of PDEs and were much better studied in the 
twentieth century. However, the analogy between odd and even-order PDEs, such as (11.11) 
and (II. 2p . is rather fruitful and will be used later on. 

The function (a nonlinear operator) g{u) in (11.21) usually corresponds to some absorption- 
reaction type phenomena and sometimes is assumed to include differential terms, such as 
D'^u, with m < 2fc + 1 (although we do not consider such cases). It is worth mention- 
ing again that, besides some special and completely integrable PDEs, general odd-order 
models such as (11.11) are less studied in the mathematical literature, than the parabolic 
even-order ones (II. 2p . 

Indeed, the most classical example of such an odd-order equation is the KdV equation: 

(1.4) Ut = Ua^xx + UUx {g{u) := UU^r), 

which was introduced by Boussinesq in 1872 together with its soliton solution [5]. The 
KdV equation models long waves in shallow water and generates a hierarchy of other 
more complicated PDEs with linear and nonlinear dispersion (dispersive) mechanisms. 
See further amazing historical aspects concerning (11.41) and related integrable PDEs in 
[23 pp. 226-229]. 

Concerning higher-order extensions, these naturally appear in classic theory of inte- 
grable PDEs from shallow water applications, including the fifth- order KdV equation, 

Ut + Uxxxxx + 30 U^Ux + 20 UxUxx + 10 UUxxx = 0. 

Let us also mention such classic examples as Lax's seventh- order KdV equation 

Ut + [35^"^ 70{u'^Uxx + uiUx)"^) + 7{2uUxxxx + 'HUxx)'^ + ^UxUxxx) + Uxxxxxx]x = 0, 

and the seventh-order Sawada-Kotara equation 

Ut [63U^ + Q3{2u^Uxx + m(^^x)^) + '^HuUxxxx + {Uxxf + UxUxxx) + Uxxxxxxlx = 0] 

see p. 234] for references and other odd-order models from compacton theory. 

Of course, PDEs such as (II. 4p and (II. ip were popular in mathematical physics last fifty 
years at least (see references below), are have been well-known in mathematical theory. 
We refer to a number of papers on local and global existence, uniqueness, smoothing, 
and various asymptotic (scattering-like) properties; see as a guide [3 [IDl [HI [121 [2Sl [33 
[SB [ini [121 HSl [S31 [SI] , where some papers on linear Schrodinger operators that exhibit 
related aspects of smoothing and evolution are also included. However, for our purposes of 
studying refined asymptotic properties of (II. ip (a refined scattering theory), these strong 
results turn out to be not quite enough. Therefore, first of all, we will need to develop 
new Hermitian spectral theory of linear rescaled odd-order operators, which will refine 
some necessary delicate asymptotic properties of PDEs under consideration. 
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In particular, as a first unavoidable step, we will need to study in detail the long- 
and short-time (on blow-up micro-scales) behaviour of solutions of the simplest linear 
dispersion equation for k = 1 (the LDE-3) 



1.5) Basic Model I : ut = Uxxx in x 



which, even in this simple standard form, will create a number of surprises that require 
difficult unusual "spectral" techniques. As a result, we then will be able to study similarity 
solutions ( "nonlinear eigenfunctions" ) of the semilinear dispersion equation 



1.6) Basic Model II : Ut = u^xx — in 1^ x where p > 1. 



These two PDEs (11.51) and (II. 6p are our main simplest canonical models to study. 

1.2. First similarity solutions: earlier history and Prandtl— Blasius results for 
odd-order boundary layer problem. For both semilinear (11.21) and linear (11.51) odd- 
order equations, we will apply the idea of self-similar solutions. In general, there are 
several techniques of solving PDEs, including using travelling wave and other group- 
invariant solutions. As customary, similarity solutions are a common way of attempting 
to understand and sometimes solve some problems. Similarity solutions are widely used 
in PDF theory, since they simplify the problem by transforming PDFs into much simpler 
ODFs, and hence this reduces the number of independent variables. Rescaled variables 
can be used in studying PDFs, since the coordinate system in which the problem is posed 
does not affect the formulation of any fundamental physical laws involved. 

The first general type of explicit solutions was traveling waves in eighteen's century 
d'Alembert's formula for the linear wave equation. The method of separation of variables 
was developed by Fourier in the study of heat conduction problems, and was later gener- 
alized and extended by Sturm and Liouville in the 1830s. Many famous mathematicians, 
such as Fuler, Poisson, Lagrange, Liouville, Sturm, Laplace, Darboux, Backlund, Lie, 
Jacobi, Boussinesq, Goursat, and others developed various techniques for obtaining ex- 
plicit solutions of a variety of linear and nonlinear models from physics and mechanics. 
Their methods included a number of particular transformations, symmetries, expansions, 
separation of variables, etc. Similarity solutions appeared in the works by Weierstrass 
around 1870, and by Bolzman around 1890. General principles for finding solutions of 
systems of ODFs and PDFs by symmetry reductions date back to the famous Lie papers 
published in the 1880s and 1890s. 

More complicated and of key importance similarity solutions next appeared in Prandtl 's 
equation, occurring in Prandtl's boundary layer theory, which was proposed in 1904, j48j . 
After that, similarity solutions of linear and nonlinear boundary- value problems became 
more common in the literature. Namely, this similarity solution in was due to Blasius 
(1908) [3] for the equation of incompressible fiuids, 

(1.7) Ipyyy + Ipx^Pyy " Ipylpxy " 'ipyt + UU^ + f/t = 0, 

where if) = il){x,y,t) is the stream function and U{x,t) is the given external far-field (at 
y = oo) velocity distribution. Remarkably, this is an odd-order PDF. 



1.3. On even-order models: the canonical heat equation. For a long period, until 
1980s, whilst many asymptotic and singularity aspects of higher odd-order models have 
not been studied in great detail, for a number of years, various even-order, both lower 
and higher order, models have been reasonably well understood. The basic linear and 
most classical case of such a PDE is the heat equation. It is one of the most important 
linear partial differential equations. A natural logic of the first part of our research of the 
LDE-3 (11. 5p can be conveniently connected with this classic area. 
Thus, the one-dimensional heat equation is the canonical PDE 

(1.8) ut = u^x in M X IR+. 
It has the classic fundamental solution 

(1.9) b{x,t) = ^F{y), with the similarity variable U — '^i 

which takes Dirac's delta as an initial function, i.e., in the sense of distributions, 6(x, 0) = 
5{x). Substituting (11.91) into the heat equation (11.81) yields a simple ODE for the rescaled 
kernel F: 

(1.10) BF = F" + \ (Fy)' = in M, J F = 1. 
This is solved explicitly to show that the solution, F, is the Gaussian: 

(1-11) ny) = ^e-4, 

which is strictly positive and does not have oscillatory components to be traced out for 
linear dispersion equations. Full spectral theory has been developed, with eigenvalues of 
the linear second-order operator B and eigenfunctions given by (see Birman-Solomjak [3l 
p. 48] for full details and spectral theory in M^) 

(1.12) a(B) = {A, = -|,/ = 0,1,2,...} and ^Pi(y) = D^Fiy) ^ Hi{y) F{y). 

Here Hi[y) denote the standard orthonormal Hermite polynomials (introduced in full 
generality by Hermite in the 1870s), which are induced by the Gaussian (II. lip . These 
Hermite polynomials in ID were earlier derived by C. Sturm in 1836 [52], where these were 
used for a classification of all types of multiple spatial zeros for solutions u{x, t) of the ID 
heat equation. This led him to the now famous Sturm Theorems on zero set] see a full 
history in [17, Ch. 1]. The Sturmian results found for the linear heat equation, have been 
the basis for other higher and nonlinear even-order and other related models, where the 
zero set analysis of solutions become indeed more involved; see [TH] and references therein. 
Whilst solutions to these models retain the basic symmetric Gaussian structure, they have 
some oscillatory tails instead of a pure simple exponential decay, which essentially change 
the meaning of Sturmian theorems. 

Thus, the discrete spectrum and Hermite polynomials in (11.120 are the cornerstone of 
the Hermitian spectral theory of the rescaled operator B in (11.101) . Since B admits a 
symmetric representation: 

^^hiipi)^ where p(y) = e^ , 
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this theory is essentially self-adjoint in the weighted space -Z^p(IR) including the A^-dimensional 
case [5]. As an unavoidable step, we will need to develop a similar theory for the LDE-3 
(11 .Sp . which is not self-adjoint and will go in different lines than the classic one. 

1.4. Very singular similarity solutions. Around the beginning of the 1980s, study of 
asymptotics of the semilinear heat equation with absorption, given by 

(1.13) Ut = Au-uP in x M, p>l (m > 0), 

led to a new class of similarity solutions called very singular solutions (VSSs). These 
are self-similar solutions, which are structurally stable in the evolution and hence, as 
t — >■ +00, attract wide classes of other more general solutions. In addition, as t — > 0, such 
a VSS concentrates at x = with infinite initial mass and this justifies the term "very 
singular" . We refer to the books [26[ HQ] for extra details and history concerning VSSs. 
Whilst VSSs were being studied early on, it wasn't until about 1985 in papers by Kamin 
and Peletier [6l [32j in which the term VSS was actually used. 
The VSSs for the semilinear heat equation (11.131) are given by 

u^{x,t) = t^p^ f(y), y = where / > solves the elliptic equation 

fA/ + iy-V/ + ^/-r = in M^ 
]f{y) has exponential decay as ?/ ^ oo. 

Existence of the VSS was established in [22] using a PDE approach, and almost simulta- 
neously, via an ODE approach, in [6]. Uniqueness of VSS was first proved later in [33], 
by using a PDE comparison method. We refer to [T^ [27j for more references and history 
of VSSs for semilinear parabolic equation. 

1.5. Back again to the KdV-type equations. It is curious that, formally, VSSs can 
be prescribed for the classic KdV equation (ll.4p . These have a standard self-similar form 

(1.14) «,(x,t) =r §/(?/), y = x/t^K 
Hence (II. 4p reduces to the following ODE: 

(1.15) /'"+ 1^/' + !/ + //' = in M. 

However, whilst we can seemingly construct a VSS for the KdV equation, we know that 
solutions of this type do not exist. The only natural solution we arrive at is the trivial 
one, / = 0. This nonexistence conclusion directly follows from the global existence results 
by Kato [551131] and Strauss [SI]. Indeed, (11.41) is invariant under refiection with 

(1.16) t \~>- T — t and x —x, 
where very singular solutions (I1.14p are then given by 

u,{x,t) = (T-t)-i/(y), y = x/{T-t)l 
Therefore existence of a nontrivial VSS would mean finite time blow-up of solutions, since 
sup^ |-u(x, t)| ~ (T — t)~i ^ +00 as t^T~. 
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In view of well-known global existence results for the KdV equation (see references above), 
this implies the non-existence of very singular solutions for fll.4l) . 

For the modified KdV equation 

(1.17) Ut = U^xx + u^u^, 

we can find VSS for p > 4, but not for p < 4. In the case p > 4, (I1.17P has self-similar 
solutions of the standard form 

u^{x,t) =t~^f{y), y = x/t^, 

where the rescaled kernel satisfies 

(1.18) f" + lyf + lf + rf' = o. 

The case p = 4 has been studied more extensively; see papers [361 M J^^^ examples. 
We thus need to study such similarity solutions for our second semilinear model (11.61) . 

1.6. Four main linear and nonlinear odd-order models to study. We will treat, in 
particular, three generalized odd-order models, which, from our "spectral-like" and refined 
asymptotic points of view, essentially have never been looked at before. In particular, we 
look for similarity solutions and using asymptotic, analytic, and some numerical meth- 
ods, we will attempt to find and justify some local and global properties of the rescaled 
solutions. We use a number of techniques, applied previously for parabolic even-order 
problems, for the odd-order ones. However, there are difficulties that arise, in particular 
due to the highly oscillatory nature of fundamental solutions and hence related linear and 
nonlinear eigenfunctions. 

Our FIRST AIM is to develop a Hermitian spectral theory of the associated rescaled 
operators, for linear dispersion odd-order equations such as (cf. (11.51) for k = 1) 

(1.19) (I) : ut = {-if+^Df^^ in R x M+ for any k > 1. 

In doing this, we gain the second goal, an understanding of its behaviour, which we 
can then use in the corresponding semilinear dispersion equation and its VSSs 

(1.20) (II) : Ut = {-l)^+^Df+\ - \u\P-\ in M X M+, where p > 1. 

In a forthcoming paper [15j, these concepts and ideas will be applied and extended to 
PDEs with nonlinear dispersion (the NDEs) of the form 

(1.21) (III) : Ut = (-l)'=+^Df +^(|u|"m) in M X M+, where n > 0. 

Concerning self-similar solutions of (ll.2ip . called "nonlinear eigenfunctions", asymptotic 
behaviour and general properties of solutions such as existence, uniqueness, shock waves, 
entropy approaches, etc., turned out to create a number of very difficult questions. First 
steps of the study of shock and rarefaction waves for NDEs including the NDE-3 such as 

Ut = {uUx)xx in X 

are performed in [201 [2^ . As a natural extension of the NDE (11.211) . we also discuss the 
VSS for the full NDE: 

(1.22) (IV) : Ut = {-lf^^Df+Wu\''u) - \u\p-\ in R x M+, where p>n + l. 
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In both cases of NDEs (11.211) and (11.221) . the key idea is to perform a "homotopy" hmit 
n ^ 0+ to arrive at the rescaled hnear operators and next to use Hermitian spectral 
theory for bifurcation-branching analysis of nonlinear eigenf unctions and VSSs. 

2. Basic linear dispersion models: LDEs and their fundamental solutions 

2.1. Basic LDE: fundamental solutions and rescaled kernels. Before looking at 
some more complicated nonlinear PDEs, it is important to understand in greater details 
how the solutions of linear PDEs behave for large and small times. This spectral theory of 
higher-odd-order linear rescaled operators will be crucial in the understanding of related 
nonlinear ones. In particular, as customary, spectral theory formed in the linear case will 
play a large role and will be used in developing understanding of bifurcations, branching, 
and asymptotic behaviour for nonlinear equations. 

Thus, we consider the corresponding linear dispersion equation (ll.lQp . which we al so 
call the LDE-(2/;;-|- 1). Whilst some lower-order cases for the odd-order linear PDE (11.191) 
such as the LDE-3 (ll.Sp are generally better understood, the higher-order cases for /c ^ 1 
are not. Indeed, it is well-known that the fundamental solution for the lowest order case 
k = 1, i.e., for (11.51) . will lead to the classic Airy function. However, we will need a more 
general calculus that is applied for arbitrary k > 1. 

Consider the self-similar fundamental solution of (I1.19P of the standard similarity form 

(2.1) b{x,t) = t^2ikiF{y), y = x/t^i^, so that b{x,0) = 6{x), 

in the sense of distributions. Substituting b{x, t) into the PDE (11.191) . we obtain the ODE 
for the rescaled kernel F, 

(2.2) BF^{-lf+^Df+^F+^^yDyF+^^F = Q in M, /F=l, 

where B denotes the first key linear rescaled operator for the LDE. Note that it is possible 
to integrate the ODE (12. 2p once, to find that F{y) solves 

(2.3) (-1)'=+^F(2'=) + = for y G M, 

which is now a linear ODE of order 2k. Further reductions of this ODE are not possible. 
For = 1, the solution of (12. 3p is the classic Airy function: 

(2.4) F" + \yF = Q =^ F{y) = Ai{y). 

Note that in our case, due to a difference in sign in the ODE, we actually have Ai{—y), 
but we shall refer to it as the Airy one. 

2.2. Asymptotic expansions of the rescaled fundamental kernel. An important 
and standard technique in trying to find the behaviour of solutions of ODEs is to use 
asymptotic analysis. This gives the limiting behaviour of solutions, in particular as y ^ 
±cx). We refer to the book by Bender and Orszag ^ for various asymptotic techniques. 
We use here a method of determining the asymptotic behaviour of the linear ODE, which 
corresponds to the classic WKBJ multi-scale analysis of ODEs, whose basic ideas go back 
to the 1920s and were reflected in a number of well-known monographs. Actually, in 
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ODE theory, asymptotics for ODEs such as (12.31) are well-known and have been classified. 
However, we will need some more refined formulae for further applications. Some of them 
for arbitrary k > 1 are not available in standard literature. Looking at the asymptotic 
behaviour of the solution of the integrated ODE (12.31) . we now write it as 

(2.5) (_l)'=+iF(2'=) = -^Fi/. 

Algebraic decay as y ^ +oo. Let us look first at the rescaled solution F{y) with oscilla- 
tory algebraic decay a.s y ^ +oo. As customary, we set, for convenience, 

(2.6) F{y) = e'^y'^ as y ^ oo. 

Now assume, as a first approximation, that s{y) is some polynomial, such that s{y) ~ ay''. 
Then the kernel F{y) and its derivatives may be given by 

F' = s'e^ F" = [s" + (s')']e^ F'" = [s'" + 3s's" + {s'fy, ... , 
p(2k) ^ ^^(2k) + . . . + ^2k - l)(s')'^-'s" + (s')2^]e^ 

where s' ~ aby''''^, s" ~ ab{b - 1) y^"^, . . . , s^^''^ ~ (b-2A:)! ^^"^^- Substituting this into 
the ODE (12. 5p . it can easily be seen that all the e'^'-^^ terms cancel. From the resulting 
equation, we look to do a dominant balance analysis, in order to determine the leading 
order behaviour of F{y). In doing so, we find a first approximation for the function s{y). 
The balance of the equation depends on the value of the parameter b and we obtain two 
different cases. 

If 6 < 0, then a.s y ^ +oo, every term on the left-hand side of the ODE is o{y). 
Therefore there is no balance in this case. If 6 > 0, then y*"" = o{y'^^^''~^^) , for any n > 1. 
Therefore balancing leading terms, we have that 

(_l)fc+la2fc^2fc^2fc(f,-l) _ __L_^, 

for all k G By first equating powers of y and then coefficients, we find our parameters 

6=^ and a = -i(2fc + l)"^(^). 
Hence, we now have the first approximation to s(y), with 

2fc+l 1 

s{y) = -2k\{^) - +c(y), s\y) = -\{^Y- +c\y), 



s"{y) = -^-^ y-'-^+ c"{y), 

2k{2k+l)^ 



2k + l , 



Here c{y) ~ o{y^^) is some function of y and the next term in the approximation of 
s{y). Since > 1, we must have that c^"^^^\y) ~ o{c^"^\y)) for any m > 1, which will be 
used in determining leading order terms. 

We attempt to find this function c{y), in order to improve the approximation of s{y). 

2k+l 

We let g{y) = —2ki\^-^^^j for convenience, to see how the terms are balanced. Hence 
balancing the leading order terms yields c'{y) ~ — ^^-^ = Integrating this, 

the next term in the expansion can be found to be c{y) ~ ^^V- The second term in 



this expansion, c{y), is called the controlling factor. It is an important term, as we will 
see that it governs the decay (or any possible growth) of solutions a.s y ^ +00. 

Whilst with the first two terms, we can have a clear idea of the leading order behaviour, 
we look to find a better approximation and now expand once again with 

2fc + l 

where d{y) = o{c{y)). Balancing leading order terms once again, we find that d'{y) ~ 
— y2 ■ Then by integrating, we find that the third term in the expansion is given by 

d{y) 

~ Note that this is the third term in the expansion and hence all following 

lower-order terms are o(l), and so, by fl2.6p . this does not affect the exponential. 
From this, we find that the required asymptotic behaviour of F{y) is given by 

2fc-l 2fc+l 

F{y) ^y-— exp{-2ki{^^) } as y ^ +00. 
Therefore, for real solutions, we have that, for some constant c G M, 

(2.7) F{y) y~^ cos [dky^ + c) as y ^ +00, where dk = 2k {^^)~ . 

One can see that (12.71) guarantees the convergence of the integral of F in (12. 2p (but not in 
the absolute sense), so that the fundamental solution b{x,t) at t = takes Dirac's delta 
in the sense of bounded measures. Note that, a priori, this is not guaranteed in general, 
since delicate ODE asymptotics are necessary to charge such a property. 

Exponential decay as y ^ —00. The same analysis can be applied for y — > — 00, by letting 
y 1-^ —y and performing the same calculations. Hence we find that the rescaled kernel 
decays exponentially fast in the opposite direction, 

2fc-l / 2fc + l N r " I 2fc+l -. 

(2.8) F{y} \y\ « cos (4|y| 2k sin 6^. + cj exp | - 4|?/| 2k j as y ^ -00, 
where dk is as in (12. 7p and 

; , , , , , I f I ^ I for even k, 

(2.9) 4 = 4|cos6,|>0, bk=\'Xl^^ „ 

U 2 +2a; for odd A:. 

We note that the coefficient a has many roots, as follows from the algebraic equation: 

r \y) ^ ^ =^ a — (-2fc+i)2fc+i 5 

which represent different solutions of the ODE. However, we only want roots such that 
there is exponential decay, rather than growth. Therefore we exclude the roots where there 
is growth, which corresponds to Re a > 0. The asymptotics here show the behaviour of 
the first roots, such that Rea < 0. The method of finding these roots and all other roots 
may be seen in Appendix A, which also explains more carefully which roots we need to 
look at. 

Thus, according to (12. 7p . as ?/ — +00, we have a purely imaginary root for all values 
of k and this gives slow decaying oscillatory behaviour. The decay of the oscillations also 
increases for larger fc, from y~^ for k = 1 towards ~ y~2 as k ^ 00. On the contrary, 
as y —00, (12. 8p shows exponential decay with oscillations for any k > 1 and these 



(2.11) \F{y)\ < 



oscillations are always exponentially fast. For the only case = 1, we have = vr, hence 
sin bk = 0. This gives the pure exponential non-oscillatory behaviour of the Airy function: 

(2.10) F(?/) ~ |y|"3e-2-3-iH4 y^_^_ 

This is the reason why the Airy function is the only odd-order linear case where there is no 
exponential infinitely oscillatory behaviour as y —>■ —oo. In a natural sense, this mimics 
the positive behaviour of the Gaussian (11. Ill) , but only partially, in a one-sided limit 
y ^ —oo. Anyway, the exponential decay in (12. 8p simplifies a part of spectral analysis of 
the operator B, for which we will need to pay the main attention to the opposite limit 
y —>■ +00, where most of technical difficulties appeared from. 

We now summarize the results of our asymptotic analysis to get a sharp uniform bound 
of the rescaled kernel F, which will eventually define the weighted functional space for B: 

Proposition 2.1. The rescaled kernel F{y) of the fundamental solution (12.11) of the LDE 
(11.191) satisfies 

' Do{l + y^)-''i^e-'^^\y\'' for y < 0, 
Do(l + 2/')~~ for y>0, 

where dk is as in {\2.f^ . Dq is a positive constant dependent on k, and 

a = ^ e{l,2) for all k > 1. 

2.3. Numerical construction of fundamental kernels. The results here were ob- 
tained by using the MatLab bvp4c solver, to look at the singular solutions of the linear 
ODE fl^ . 

Obviously, if F is a solution of (12. 3p . then cF is also a solution, for all c G M. Due to 
such "non-uniqueness" (and, to some extent, an "instability") of the solutions, the zero 
solution is likely to be found using numerical methods. Therefore in order to ensure that 
a non-zero solution is obtained, we set the (normalisation) constraint 

(2.12) max|F| = l, 

which is attained at some point y = a. We then solved the ODE for different right and 
left solutions at this maximum point a, using the bvp4c solver. Since we have fixed a as 
a maximum, the first derivative is also zero at this point. The other boundary condition 
placed was to ensure that F{ii) = 0, at an end point that is sufficiently removed from 
a. In order for a match of right and left solutions at the point a, we needed to have the 
second derivative of the solution to be the same here, in order for F{y) to be continuous 
at y = a. So the value of d was moved in order to match the second derivative for the 
right and left solutions. See Figure [U which was obtained by this shooting method. 

Indeed, for /c = 1, as a simple alternative, one can use "shooting from the left" in the 
second-order ODE (12.41) by using the exponential decay asymptotics (I2.10p for y ^ —1 to 
get a proper behaviour (12.70 for y ^ 1, since no other are available in the version (12. 4p . 
Such a simpler one-sided shooting is shown in Figure [2], which is rather similar to Figure 
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([T])(b). What is most important is that the sizes of the tails in both Figures coincide. 
However, such an approach is hardly applied in higher-order cases k > 2, where extra 
solutions selecting procedures such as mult i- dimensional shooting are unavoidable. 

A similar method of two-sided shooting was applied to the fifth order {k = 2) equation, 
where values of the left-hand solution were used as boundary conditions for the right 
solution, and the correct value for d was found by matching the fourth derivative. See 
Figure [31 Note that, in view of fl2.12p . the plots do not show the fundamental solutions 
of the ODE, such that J F = 1, but rescaled profiles, for which the normalization (12.121) 
is satisfied. 

As one can see from comparing Figures [T] and [31 the oscillations are smaller as k 
increases. This follows from the asymptotic analysis done in Section fI72[ see fl2.7l) . Decay 
of the algebraic envelope is very slow as y ^ +oo also, especially for the Airy function 
when k = 1, and is faster for A; = 2 in Figure [31 

For convenience we denote these rescaled kernels as higher-order Airy functions by 

(2.13) F{y) = Ai2k+iiy) for A: = 1, 2, . . . , so that Ai = Aig. 




-20 -15 -10 



10 15 20 




-20 20 40 60 



100 120 140 160 180 200 



(a) Oscillations of F{y) for y e (0, 20). 



(b) Oscillatory "tail" of F{y) for y > 1. 



Figure 1. The rescaled kernel F{y) = Ai(y) of the fundamental solution (|2.1 
of (jl.Sp (k = 1) obtained by two-sided shooting. 



3. Explicit semigroup representation: first spectral properties 

Here, we begin to develop first aspects of a Hermitian spectral theory, which our further 
analysis will depend crucially on. 

3.1. Operator B: formal series. Let u{x,t) be the solution of the Cauchy problem for 
the LDE-(2/c -|- 1) (ll.lQp . with bounded measurable initial data u{x,0) = uo{x) from a 
weighted space to be specified. The solution is then represented by the convolution of 
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Rescaled Airy kernel F(y), k=1 




Figure 2. The rescaled kernel F{y) = k\{y) of the ODE (p^i]) {k = 1) obtained 
by one-sided shooting via the bundle (|2.10p . 
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Figure 3. The rescaled kernel F{y) = A\^{y) of the fundamental solution (j2.1 
of the LDE~5 Ut — ~Uxxxxx — 2) obtained by two-sided shooting. 



initial data with the fundamental solution given by a Poisson-type integral: 
(3.1) u{x, t) = b{t) *uo = / - z)r^)uo{z) dz. 



Let us next introduce further rescaling corresponding to the variables of the fundamen- 
tal solution fl2.ip . 

(3.2) u{x,t) = t~^i^w{y,T), where y = x/t^i^, r = Int : R+ M, 
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where we now have scahng with respect to time as well. The rescaled solution w{y, r) 
then satisfies the evolution equation 



(3.3) Wr = Bw, 

where B is the linear operator for the rescaled kernel described in (12. 2p . More precisely, 
B is still a formal linear differential expression to be equipped with proper "boundary 
conditions" to get its actual domain. 

Here, w{y,T) satisfies the Cauchy problem for (13. 3p in M x R+, with initial data at 
r = 0, i.e., ai t = 1 and not t = 0. So now the initial data are given by 

(3.4) wo{y) = u{y,l) = b{l)*uo = F*uo. 

Hence, the linear operator ^ — B is the rescaled version of the original linear disper- 
sion operator ^ + (— 1)^D^.^+^. Therefore, the corresponding semigroup e^"^ admits an 
explicit integral representation. This helps to establish some properties of B, including 
key spectral ones, and describe other evolution features of the linear fiow. 

Thus, rescaling convolution (13.11) gives the explicit representation of the semigroup 

(3.5) w{y,T) = J F[y — ze^'^)uo{z) dz = e^'^u{y,l) for r > 0. 

For any |/ e R, Taylor's power series for the analytic kernel F can be used to expand the 
convolution and to obtain 

(3.6) F{y-ze~^^)=j:e'^^i^DP^F{y)z^^Ee-^^^My>^- 

Here, for the first time, we reveal the eigenvalues and eigenf unctions of the operator B as 
(3-7) = - J^ and ^p{y) = ^-^^^^D^F{y) for any (multiindex) (3, \(3\ > 0. 

Remark: on extensions to M^. We note that whilst \(3\ = j3 here, this is not the case in 
the multi-dimensional M^, where (3 stands for a multiindex. However, we often still keep 
the notation of to show that the theory here can be extended to multi-dimensional 
spaces. We bear in mind that Hermitian spectral theory can be properly developed for 
the typical linear dispersion operator in M^, 

(3-8) B = i-ir'i-A^ + wriy'^y + wri^^ 

which appears after the scaling (similar to ( 13. 2p ) of the following LDE: 
(3.9) ut = {-l)'£^A>'u, 

where analogies of the operator on the right-hand side occur in completely integrable 
equation theory (cf. the Kadomtsev-Petviashvili equation for N = 2 and others from 
higher-order hierarchies). The multi-index is then /3 = {f3i, . . . , (3^), with the length 
\j3\ = Pi + . . . + Pn- With this notation, some of our basic results can be directly 
translated to operators such as (13. 8p in R^; see [131 [H] for necessary technical details. 

13 



Thus, without loss of generahty, we continue to develop a ID theory. Looking back at 
the expansion of the convolution fl3.5l) via fl3.6l) . the solution of (13. 3p is given by the series 
(understood formally still): 

(3.10) w{y, r) = ^ e'^^ Mp{uo)^M- 

It then follows that \p = —i^^^ ^^^d il^piy) are the eigenvalues and eigenfunctions of B as 
in ([32D and 

(3.11) Mp{u^) = ^J zPuo{z)dz, 

R 

then turns out to be corresponding moments of the initial data Wq, i.e., "scalar products" 
of with some "adjoint polynomials" (eigenfunctions of the "adjoint" operator B*) to 
be detected shortly together with a proper metric involved. 

3.2. B: convergence of series. We next need to check suitable metrics of convergence 
of formal series (13.101) . Those questions were addressed in [13] and in the most general 
case in [18] in the purely "elliptic" ("parabolic") case, when the rescaled kernel F{y) 
always have exponential decay as ?/ — oo. In view of (12. 7p . this is not currently the 
case as y +oo. We concentrate on this case and refer to the analysis in [TH § 5] of 
convergence of the integrals as y —oo. 

First of all, as in [THl ITS] , we have that the ODE (12.21) for F implies the following 
estimate of the eigenfunctions: 

(3.12) |^^(^)|<-^(l+^)^-2|^«f for y>0, where «fc=^4sin6fc. 







Here, the main growing factor ~ |/2fc as ?/ — >■ +oo, as well as the multiplier a^, come from 
differentiating f3 times inside the cos function in (12.71) . 

We next claim that the series in (13.101) converges uniformly on compact subsets for data 
from the weighted space (e.g., with fast exponential decay at infinity): 

(3.13) uoeLjiR), where p(y) = e'^l^'l", a = 

and a > is a sufficiently small constant. As in the parabolic case [13], in view of the 
exponential decay (12.81) . for ?/ < the condition a G {0,2dk) is necessary, but this case 
will not be treated here as being standard. 

The proof of uniform convergence on compact subsets uses typical Holder estimates 
such as (for y > again) 



(3.14) J\zf\uo\^Jj.\z\^Vp\uo\ < V/pKP J II 



1 2/3 



and, using standard properties of the F-function and Stirling's formula, we estimate the 
last integral as follows: 



(3-15) /4W^~r(ii)~(ii),^(_|M_)^ 
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Overall, this gives the following majorizing series for fl3.10p for y G (0, L): 



(3-16) E(,);^V(^)'^^' 

which obviously converges for arbitrary fixed L > 0. 

We next show that the series (I3.10p converges in the mean in the metric of the weighted 
space 

{pilyr for w < —1 
e y for y > 1, 

cf. f l3.13p . Using the same estimates (13 .120 and (13.140 . we then derive the following more 
sensitive geometric majorizing series: 

(3-18) E(,)7r^, where 7. = ^^(^Itt)''"' < 1' 

so the series converges. We do not include in 7^ the extra multipliers ~ \dksmhk\ < 1 
shown in (I3.12p , since the first one < 1 already suffices for the convergence of (I3.18P . 

3.3. B: alternative representation of semigroup and series, generalized Hermite 
polynomials. Let us discuss an equivalent explicit representation of the semigroup for 
B, which, though being more difficult, will more clearly determine the eigenfunctions of 
the adjoint operator B* to be introduced and studied next by a simpler direct approach. 
We perform another rescaling to exclude the relation (13. 4p . in order to find the correct 
semigroup, corresponding to the initial data at t = 0, i.e., Wq = Uq\ 

u={l + t)-^w, y = x{l + t)~^, r = ln(l + t) : M+ ^ M+. 
Then rescaling the convolution gives 

(3.19) w{y, t) = e^^Mo = (1 - e~^)~2FrT J F[{y - ze~^^){l - e^'^y^^)wo{z) dz. 

Once again we can look to find explicit representations for the eigenfunctions, eigenvalues, 
and adjoint eigenfunctions, given in the dual products {uoyip^), in the standard metric of 
(actually, a more delicate indefinite metric will be needed; see below). It will be shown 
later that we can actually determine the adjoint eigenfunctions, {^^}, using a much easier 
method. 

Meantime, let us perform some convenient manipulations. Looking at our rescaled 
equation (13.190 . by Taylor's expansion we have 

F{{y - ze-^){l - e"^)"^) = J2(f^) - e-^)"^)(e^ - ly^zf" 

and F{yil - e-^"^) = Eh ^ (DyF^my^il - e-y^. 
Then, using these expansions, our solution is given by 

^3 ^) = (1 - e-^-^ {D-rmD^iynie^ - 1)"^ 

x(l - e~'^)~2^ J^z^wo{z) dz. 
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Rearranging this, we have 

(3.21) wiy, r) = .) e-^(l - e-^""^ ^-0' {D'^FmD^iyn I^z^wo{z) dz. 

Further expansion in the exponential term containing (1 — e^^) ^fc+i in terms of e~^ 
yields an alternative more complete eigenfunction expansion of the semigroup with a full 
representation of all the eigenfunctions and the generalized Hermite polynomials {V'^} in 
the moments M^. Moreover, being re- written in the form 

(3.22) w{y, t) = e~^^p{y) K, V-;)*, 

this expansion, besides the adjoint polynomials $* = also defines the corresponding 

indefinite metric (■, to be carefully introduced and explained in Section [5l where those 
polynomials will be obtained via a simpler direct approach. 

3.4. Classification of large-time asymptotics. As a by-product, for our LDEs, we 
have described the large-time asymptotic behaviour for the problem 

(3.23) ut = i-lf+^Df+^u in M X M+, where Uq E Lj, 
and the weight p is as in fl3.13p . 

Theorem 3.1. For the problem (13.231) . the eigenfunction expansion (I3.10p implies that, 
for any uq G L^^(R), uq ^ 0, there exists a finite I such that, as t +oo, 

u{x,t) =t'^^^^[cf3^f3{x/t^) +o{l)], where ^13 = -^ 

and I = \(3\ is the first eigenvalue index, with the corresponding moment cp = M^(mo) 7^ 0. 

The following uniqueness conclusion is straightforward and we keep this as a simple 
illustration for further results: 

Corollary 3.1. Assume that, for any constant K > 0, the solution of fl3.23p satisfies 

sup^ \u{x,t)\ = o{t'~^) as t ^ 00. 

Then u{x, t) = 0. 

Such results belong to classic Carleman-Agmon-type estimates in operator theory: if 
a solution of a linear equation, under proper conditions on operators involved, decays 
super-exponentially fast (in terms of r = Int) as t +00, then it is trivial. For elliptic 
equations P{x, D)u = 0, this has the natural counterpart on strong unique continuation 
property saying that nontrivial solutions cannot have zeros of infinite order; a result first 
proved by Carleman in 1939 for P = -A + V, V e L^^, in [8J; see [SDl [HI |55j for 
further references and modern extensions. 

Thus, for convenience, we fix again some conclusions achieved above: 

• There exists point spectrum {X/b = — 2Ah^' — 0} of non-symmetric operator B. 

• We have no integral terms in the expansion, hence the spectrum is expected to be 
discrete in a proper functional setting. 
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• The set of eigenfunctions {ip/3} seems to be complete and closed in a certain 
weighted L^-space. 

• Traces of the polynomials, which give rise to the adjoint eigenfunctions, {V'^}. 

3.5. Semigroup of the adjoint operator B* and its eigenfunction expansion. We 

now find the explicit representation of the semigroup e^*"^, where B* is obtained from the 
LDE in (13.231) by using other blow-up rescaling. Let us introduce the rescaled variables 

(3.24) u{x,t) = w{y,r), y = x/{l-t)^^, r = - In (1 - t) : (0, 1) ^ R+, 
where, for convenience, the blow-up time is T = 1. Then w now solves the problem 

Wr = B*iy for r > 0, 
with data w{y, 0) = Wo{y) = Uo{y). Here the "adjoint" operator B* is given by 

(3.25) B* = {-ir-^^Df^^-^^yDy. 
By rescaling the convolution (13. ip . we have 

w{y,T) = (1 -e"^)"^ J^F{{ye^^ - - e-''y^)wo{z) dz. 
Using Taylor's expansion yields 

F((?/e"2fcTT - z){l - e"^)~2fcTT) 
=E(/3)^'^(1 - e~^)-^D'^F{ye'^{l - e-^)-^)zf^ 
and expanding in y leads to 

F{ye-^{1 - e--)-^) = ^D"^ F{0)y''{e^ - ly^. 
So the solution is represented as 

xD^FiO)j^y'^-^J^zPwoiz) dz. 

Again, similar to the expansion (13.221) for the operator B, (13.261) can be viewed as an 
eigenfunction expansion of the solution that can reveal many key spectral properties of 
B*. In particular, we obtain the same real spectrum and eigenfunctions: 

(3.27) B* : A/3 = —21^:1 and Hermite-type polynomial eigenfunctions ip'^iy)- 

However, further refining of this expansion will lead to more complicated formulae, which 
are not that effective and useful. Therefore, we return later on to polynomial eigenfunc- 
tions of B* on the basis of a simpler direct approach. 
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3.6. Classification of multiple zeros: micro-scale blow-up asymptotics. Due to 

the blow-up character of the scaling (13.241) . the polynomial eigenfunction set $* = {V'^} 
of the adjoint operator B* is able to give insight into the structure of multiple spatial 
zeros of solutions u{x,t) of the corresponding LDEs. This blow-up asymptotic theory is 
a natural counterpart of the large-time one in Section 13.41 induced by eigenfunctions $ of 
B. The key principles of such a classification for LDEs can be found in fiSi § 9], where 
other blow-up application of such a spectral analysis can be found. 

4. HERMITIAN spectral THEORY: OPERATOR B 

We now start more systematically to develop necessary spectral theory for the operator 
pair {B, B*} introduced above. In some aspects, this theory repeats standard steps of self- 
adjoint theory for the Hermite classic operator fll.lOp . which since the nineteenth century, 
is associated with the names of Sturm, Hermite, and other famous mathematicians. 

4.1. Bounded operator B. According to (13.171) and (I3.13p . we calculate the spectrum 
of the linear operator B in the weighted space L^(R), with an exponential weight given 
in (I3.17p . where, as usual, we assume p{y) to be sufficiently smooth in the complement 
interval [1, —1]. Here we have that p{y) > and a G (0, 2dk) (this is necessary for y <^ —1 
only) is a sufficiently small positive constant, where is defined as before, in (12.70 . (12. 9p . 

We introduce a Hilbert space of functions iJ^'^"^^(M) with the inner product and the 
induced norms: 

2k+l 2k+l 

{v,w)2k+i,p = J p{y) E Dlv{y)Dlw{y)dy, \\v\\lk+i,p = J p{y) E \Dlv{y)\^dy. 

It turns out that it suffices to restrict to the real case. Then H'^^'^^ C Iv^(R) C L^(R). 
The first conclusion is standard as in 



Lemma 4.1. B is a bounded linear operator from Hp''^^(R) to L^( 
Proof. Consider the differential expression (12.20 . 



Bv = (-l)'=+li;(2fc+l) + 1 + 



V. 



2k+l y ' 2k+l 

For B to be bounded, it is necessary to look at the second term with the unbounded 
coefficient y. In order to do this, we want to show that 

Jpiyv'Y < C/p(t;(2A^+i))2d?/, 

for some constant C > 0. To show this, we look at the non- negative integral 

< / p{v' + y^vf dy = J p{{v'y + y'^^v'^ + 2y^v'v) dy, 

where 7 > is some unknown exponent. Integrating by parts (allowed for functions in 
H'f'^^iM)) in the last term yields 

2j{py^)v'vdy = j{py^){v^ydy = - j v\py^y dy. 
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Then for the exponential weight, we have to have that 

p{y) = e-'^y" =^ {py^y ~ e-"^'^y^+"-\ 
for y 3> 1 (and also for y <^ —1 that is similar, where we replace y ^ \y\). Hence 

/ p{v'y dy + J py^'<v^ dy + Cij py^+^^^h^ dy > 0. 
By equating powers of y, this yields 7 = a — 1. Substituting this 7, we have that 

J py^i'^-^^ dy < C2 J p{v'y dy. 
In particular, for = 1, by the Hardy-type inequality, the following holds: 
J py^(»~^){v'ydy < Cs J py^^'^-^^vy dy < Cijp{v'ydy. 
However, we want that 

Jpy\Frdy<C,JpiF'ydy. 

Hence we have 4(a — 1) = 2, which gives a = §• 
For the general case, we have 2k iterations and so 

4k{a-l) = 2 =^ a = □ 

4.2. Completeness and compact resolvent. The following conclusions are also pretty 
general jlSl ITS] : 



Lemma 4.2. (i) The set of eigenf unctions $ = {^Jp} given in ^3. ?D is complete in 
and in i^^(M). 

(ii) The resolvent (B — A/)~^ is a compact operator in L?(]R). 



Proof, (i) In fact, due to the differential structure of the eigenf unctions in (13.71) . this 
is a standard result in functional analysis; see analytic function theory approach in 
Kolmogorov-Fomin [381 P- 431]. Let us show that the system of eigenfunctions {D^F} is 
complete in L?(M); the case of L^(M) is simpler. By the Riesz-Fischer theorem, we have 
to show that, given a function G G L'j^iM), then 

/ Df^F{x)G{x) dx = for any 

implies that G = 0. Let F = T{F) and G!F{G) be the Fourier transforms of F and G. 
Then substituting F = J^'^{F), G = J^-\G), and using that D^J^-^{F) ~ ^^J^-\F), 
after integration by parts via JF^^(G')(^) ~ JF(G')(— ^), we obtain that 

(4.1) /e^^(OG'(-Ode = for any /?. 

Let us calculate F. Applying the Fourier transform to the equation (12.21) yields 

(4.2) i\^\2k^P + _l_^DF = ^ F(0 = e-'l«l"«. 
Substituting (g^D into (O) yields 

(4.3) /e'^e-'l«l''«G'(-Ode = for any /3. 
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Finally, consider the function 

M{z) = /e-'l«l''«G(-Oe'^Me, 

which is entirely analytic in C, [37j. So fl4.3p means that D^M{0) = for any (3. Therefore, 
M{z) = 0. Hence Q{C,) = almost everywhere and G = 0. 

(ii) The proof follows that of the 2mth-order case, set out in [TSj § 2.3]. In fact, replacing 
m I— i> k + ^, will yield the same result. This also directly follows from the compact 
embedding of the corresponding spaces H'^^^^ C L^; see Maz'ya's classic monograph on 
Sobolev spaces [HI p. 40]. □ 

4.3. Space of closure L^(R). We next apply another direct approach to constructing 
the necessary domain of B; cf. that in [13]. Such a rather "artificial" construction, though 
not being completely unavoidable for such odd-order operator B (a more traditional way 
of doing this along the lines in [T3| [T8] is still available at least partially), will essentially 
simplify in Section [5] the way of introducing the set of linear functionals $* = {V'^} from 
the adjoint space being generalized Hermite polynomials as eigenfunctions of the adjoint 
operator B*. 

Namely, as in [TH § 5], given the complete set of eigenfunctions $ = {ipp} we introduce 
the space of ^-closure denoted by i^p(IR) C -^^^(M) as follows: 

(4.4) Lj(M) : closure of finite sums 'Ylip=o^l3'4^P metric of L^(M). 
The notation the "space of closure" is then justified naturally: 

(4.5) any v G I'p(IR) has the unique eigenfunction expansion via $: v = X]{/3>o) 

i.e., the set $ is closed therein. Note that, for the whole space i^p(IR), the property (14.50 
is not easy to check and, most plausibly, it fails. A functional characterization of such 
a space is not that straightforward, and some estimates of the sequences {c/j}, for which 
the series in (14.51) converges, can be derived as in [rF, § 5], so we will not treat those 
questions here. In Appendix A, we present an alternative explanation of some of key 
"radiation condition" properties for functions from the space of <l>-closure. Appendix B 
shows spectral theory for a majorizing integral operator for (13. ip . which is order-preserving 
so some comparison features are available. 

We now complete our spectral theory of B as follows: 

Lemma 4.3. The operator B m f TOj) with the domain V{B) = Ll{R) n H'^''+'^{R) has 
the point spectrum comprising real eigenvalues only: 

(4.6) a(B) = {A^ = -2^,/5 = 0,l,2,...}. 

Eigenvalues are simple with eigenfunctions [3.1). 

Proof. Eigenvalues and eigenfunctions are found by applying to (12. 2p : 

By the construction of the space of closure (14.40 . this completes the proof. □ 
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Note also that as follows from the asymptotic expansion of (13.101) and (13.201) as r — >■ 
oo, no other eigenfunctions exist (at least for such data). This again confirms that all 
eigenvalues are real and are given in (14. 6p . 

A sharper using Stirling's series shows that (see estimates in [18^ § 5.1] that are applied 
to the present case) 

(4.7) = ^c^^^ e L2 ^ = o(|/5|l'^l('^+")) withan£>0, 
for "almost all" ^ 1. A sharper and almost optimal estimate reads: 

(4.8) Ecp^P^eLl if cp = 0{6\^\\f3\\^\n, where u = '-^ = 

with a sufficiently small constant 6 > 0. 

By H"^^^^ C I/^, we denote the dense linear subspace obtained as the closure in the 
norm of H^^^^ of the subset of eigenfunction expansions with coefficients satisfying (14.71) 
or (14.81) . Hp'''^^ with the scalar product of H"^^^^ becomes a Hilbert space and can be 
considered as the domain of B in if^'^^^. There holds 

(4.9) iff +1 C p 

4.4. Little Hilbert spaces. We will need a subspace of introduced as a little Hilbert 
space of functions f = ^ cpipp G with coefficients satisfying 

(4.10) E|c/5p<oo, 

where the scalar product and the induced norm are given by 

(4.11) {v,w)Q = Y,cpO'P ^oi w = Y,api)p ePp, and \\v\\l = {v,v)o. 

Obviously, is isomorphic to the Hilbert space P of sequences {cp} with the same inner 
product, and hence 

(4.12) $ is orthonormal in Z^. 
It follows from ( 1T8i) that 

(4.13) lldLl =^ llr\Ll = ll. 

We next define a little Sobolev space h'^J'^^ of functions v E Ip such that Bu G Ip, i.e., 

^|A/3C/3p < OO. 

The scalar product and the induced norm in /i^'^'^^ are 

(4.14) {v,w)i = {v,w)o + {Bv,Bw)o and \\v\\l = {v,v)i = J^i'^ + \Xp\^)\cp\\ 

This norm is equivalent to the graph norm induced by the positive operator (— B + al) 
with a > 0. Then h"^^^^ is the domain of B in I p. We also have a Sobolev embedding 
theorem, 

(4.15) hf+^ C ll compactly, 
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which follows from the well-known criterion of compactness in P: a T C P is compact iff 

(4.16) V £ > 3 integer K = K{e) > such that V {cfs} G T ^ T.\(3\>k < ^; 
see, e.g., |13]. 

Since the orthonormality of $ is known to be of importance in operator theory and 
applications, in some linear and nonlinear problems dealing with operators like B, the 
little space can play a special role in comparison with the big one L^. 

It follows from (14.111) that B is self-adjoint in with the domain h'^J''^^, 

(4.17) (Bv, w)o = {v, Bw)o for all v,w e hf^\ 

Notice that this a posteriori conclusion in a special functional setting is obtained after 
establishing all the necessary spectral properties of the operator. 

5. Spectrum and polynomial eigenfunctions of the adjoint operator B* 

We now look to explicitly describe the eigenfunctions of the "adjoint" operator (I3.25p . 
We recall that B* was derived in Section via the blow-up scaling (13.241) . and yet has 
nothing to do with the standard adjoint differential form (B)* in the metric of L^, which, 
though is easily derived, has no further applications. 

5.1. Indefinite metric. Thus, before we look at the operator B* however, we first make 
the following easy observation: 

Proposition 5.1. B* is not adjoint to B in the standard metric o/L^(]R). 
Proof. Let v,w & C^(]R), then integration by parts yields 
(B^, w) ^ /((-l)'=+i^;(2^-+i) + ^ {yvy)wdy 

R 

R 

= ... = /^((-l)V2'=+i)-2^y^') dy 

R 

= {v, (BYw), where (B)* = (-1)^=^^ +i --^yDy^ B*. □ 

Thus, in order to get the correct adjoint operator B*, it is necessary to use another 
metric. The scalar product of this indefinite metric of the space is given by 

(5.1) (w, w)^ = J v{y)w{-y) dy. 

R 

Since B and B* have real point spectrum (see (I3.27p . a full theory to be developed), we 
may omit the complex conjugate. 

Proposition 5.2. B* is adjoint to B in the indefinite metric o/L^(]R), with the indefinite 
scalar product (15.11) . 
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Proof. For our operator, taking v,w & C^{M.) 

R 

= J((_l)fc+V2fc)(^) + _l_^^)«;'(_^)dy 

R 

= ... = Jv{y) _ _L_ (-y)^'(-y)) dy = {v, B*w),. □ 

R 

Thus, B* is adjoint to B in the given indefinite metric, which we write down again as 

(5.2) {v,w)^ = f v{y)w{—y)dy = {v, Jw), for any v G L^, Jw G Lf^^. 

Here, the canonical symmetry operator Jw{y) = w{—y) is bounded, self-adjoint, and 
unitary (it is the Gramm operator of this metric). Moreover, one can see that, for the 
given anisotropic weight p{y) in (13.171) . 

Jw G Ll/p, if w G LJ, 

so that the adjoint space in the indefinite metric (15. 2p is 

(5.3) (Lj): = Ll 

Let us next introduce another natural weighted space, with a symmetric exponentially 
decaying weight, 

(5.4) L2.(R), where p* {y) = = e-\y\" in M, so (L?)^^ = L^.. 

Note that the standard "adjoint relation" of the weights p*{y) = exists for y <^ —1 
only and is wrong for y ^ 1. Indeed, this refiects a strong anisotropic behaviour of the 
eigenf unctions 'ippiy) of B as ?/ ^ ±oo. 

The set of even functions Ej^ = {v{—y) = v{y)} is a positive lineal (a linear manifold) 
of the metric (15.21) , 

(f , v)^ > for V e E+ C Lj, v ^ 0, 

and odd functions = {v{—y) = —v{y)} C give the corresponding negative lineal. 
Therefore, with this metric is decomposable: 

v = v+ + v. = "(^)+;(-^) + -iy)--(-y) ^ where v± e E± ^ Lj = E+ ® E_, 

where, in addition, E+A^E^ in the metric (15.21) . The corresponding positive majorizing 
metric is given by 

\{v,v)^\ < [v,v]^ = {v+,v+)^ - {v-,v^)^, 

etc. This case of the decomposable space with an indefinite metric having a straight- 
forward majorizing one is treated as rather trivial; see Azizov-Iokhvidov [T] for linear 
operators theory in spaces with indefinite metric. Metric (15. 2p is widely used therein; see 
[II p. 13, 17, 23, 114]. Then the domain of B* is defined as H^^,, etc. 

Historical Remark: As we mentioned, basic results of linear operator theory in spaces 
with indefinite metrics can be found in Azizov and lokhvidov [1]. It was not until about 
1944 that L.S. Pontryagin published the article on "Hermitian operators in spaces with 
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indefinite metric" [U]. A new area of operator tlieory liad been formed from Pontryagin's 
studies, wliicli, during tfie time of the WWII, were originated and associated with some 
missile-type mihtary research, [H]. This work set by Pontryagin was continued from 1949 
and in the 1950s by M.G. Krein [39] and I.S. lokhvidov [29j. 

5.2. Discrete spectrum and polynomial eigenfunctions of B*. We consider the 
spectrum of the linear adjoint operator B* in the weighted space L^*(M), with exponen- 
tially decaying weight (15. 4p . Here we have that p*{y) > and a G (0, 2dk) (see (12. 9p ) is a 
sufficiently small constant. The proof of the following results does not differ from that of 
Lemma 14. 1[ 

Lemma 5.1. B* is a bounded linear operator from iJ^^^^(]R) to L^, (M). 

As we have mentioned before, we now, in an easier and direct manner, construct poly- 
nomial eigenfunctions of B*. 

Lemma 5.2. (i) The operator 25\i in L^, has the same point spectrum o"(B*) = a"(B) 
as in ( |^.6| ) and the eigenfunctions ip* (y) are generalized Hermite polynomials given by 

(5.5) rf^iy) = ^ y^+i-l)"^' E j^D^^'^'^^ yf^ . 

(ii) The eigenfunction subset $* = {ip^} is complete in and in L^, . 

Proof, (i) The construction is straightforward using the fact that the operator yDy is 
homogeneous. Consider the eigenfunction equation 

(5.6) B>; = xrp. 

Assume that ip^{y) is a polynomial of degree Z > 0. Then we look for a solution of (15.61) 
in the form 

where each Pj{y) is a homogeneous polynomial and s = [afc+iJ- From the eigenvalue 
problem (15. 6p . we can work out all terms of the polynomial, for a given degree /. Then 

(5.7) (-1)^+1^2^+1^; -^^D^^; = A;3V^; =^ A;3 = -^, / = 0,1,2..., 
and we can define all other polynomials Pj{y) by 



(5.8) PM = ^^^D(''^''^^Po{y), j = l 



(fc+i) 



s. 



Fixing the leading term Poiy) = y^ calculating all the terms via (15. 8p yields (15. 5p . 
(ii) The proof of completeness is much simpler than that in Lemma l4.2( ii). In fact, it 
is a common fact in functional analysis that polynomials are complete in any reasonable 
weighted spaces; see the famous text-book [3Sl p. 431]. □ 

As in (14. 4p . we next define the space of $*-closure L^*(]R) with the corresponding 
property (14. 5p . Since, obviously, the resolvent now acts as follows: 

(5.9) for^ = Ec/3^;, (B*-A/)-it; = Ei^^^, 



'^13- 
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i.e., the convergence of the series is improved, it is not difficult to conclude that the 
resolvent is a compact operator in L^*, though we are not going to use such facts. By 
construction, fl4.6p is the full spectrum of B* in L^, , with the domain 2^(6*) = L^*(]R) fl 

We need a few more properties. Firstly, we note that, as in [HI § 5.2], given an 
arbitrarily small e > 0, there holds 

(5.10) = o(i/?ri/5i(.+.)) =^ Y^cpr^^Li.. 

In addition, if t> = '^Cisip*^ G L^. , then for arbitrarily small fixed e > 0, 

(5.11) = od/JI-I^K'^-)) for "almost all" » 1 (^=^ = 2^)- 

In view of the fast decay (15.101) of the expansion coefficients, similar to we introduce 
the adjoint little Hilbert space of eigenfunction expansions v = Y2 cpipp ^ -^p* with the 
scalar product (■, ■)o* and the norm || ■ ||o* defined as in (14.111) . Then the domain of B* 
in is the corresponding little Sobolev space /z^*^^ compactly embedded into Z^*. Here, 
(-, ■)i^, and II ■ ||i=K denote the scalar product and the induced norm. Then B* is self-adjoint 
in , and L^, , Hpt~^^ are dense subspaces of 1^, . 

5.3. Bi-orthonormality of the bases $ and $*. This is a principal issue for several 
applications of eigenfunction expansions via eigenfunction sets $ and the adjoint one $* 
to be developed later on. First of all, we note that the generalized Hermite polynomials 
'ippiy) G initially appeared in the semigroup expansions (I3.10p . (13.111) and (I3.22p . as 
more or less standard linear functionals defined in the weighted space i^^(M) as given in 
(13.131) . Therefore, by Hilbert space theory, these polynomials belong to the true adjoint 
space: for any /3, 

(5.12) rpeLf{R) = Ll,{m), where p* = j (cf. (El). 

Moreover, the whole set $* of such polynomials is complete in L^, , and is closed in the 

corresponding space of $*-closure L^*. 

However, this is not enough for further applications. Namely, let us return to the space 
of ^-closure i^p(M) defined for B in Section 14. 3[ According to the construction therein, 
we have defined linear functionals from the adjoint space L^*. Thus, as we have seen 
earlier from the eigenfunction expansion of the semigroup (13 101) and (13.211) (for r = 0), 
those linear functionals are well defined for data uq G L'j with the standard (not in the 
v.p. sense) definition of the integral; cf. (13.110 . Therefore, we then need to extend those 
functionals to the whole space of ^-closure. 

Before doing that, let us get rid of the indefinite metric (■,■)*, applying of which could 
rise some natural questions, and return to the standard (dual) L^-one. Actually, we then 
need slightly revise our definition of eigenfunctions {'0^(?/)} of B*. Looking back at the 
expansion of the convolution, where the rescaled solution is given by (I3.10p . we note the 
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definition of the moments of the initial data, which give rise to the adjoint eigenfunctions. 
These moments f l3.1ip . in the L^-metric, are assumed to be given by 

Mf3{uo) = {wo{z),ip*p{z)) = ^ J zf^uoiz) dz for any p. 

However, we must have the representation in the indefinite metric that actually yields 

Mpiuo) = {wo{z),rp{z)), = f{~zyuo{z) dz=^-^ fz^Uoiz) dz. 

Therefore, in the expansion there must be an extra multiplier (—1)'^ for the dual prod- 
ucts. So, for purposes of convenience, instead of (15.51) . we now have to take the adjoint 
eigenfunctions in the following form (i.e., the extra multiplier (— l)'^' is included): 

l^-l 

(5.13) rpiy) = ^-^\y'+i-iY'^'^ E j^D(''^')^y^ 

L j=i 

Thus, as the next step, according to our construction, we first note that, by standard 
properties of Hilbert spaces, the space adjoint to L'j^(M.) is 

(5.14) Ll*{R)=L%R), where p = ^ ^ P* =^ rp^LfiR). 

Similarly, in the metric (15. 2p . for which (15. 3p holds, there is no inclusion of polynomials 
either. In other words, unlike the more standard "parabolic" case [I3l [18] (where always 
p* = i), the polynomials ip'^ do not belong to the adjoint space. However, we deal with 

the space of closure C L^, with the adjoint one L^* being wider than the pure and 
standard L^*. Finally, completing this discussion of such intriguing aspects of functional 
spaces involved, which will not play a role later on, we mention that one needs to extend 
and define such "extended" linear functionals {ip*^} from the adjoint space L^* to the space 

L^, , which are defined in the whole L^. On one hand, this looks related to a standard 
procedure of extension of uniformly convex functional in linear normed spaces by the 
classic Hahn-Banach theoren£|, [38j . 

On the other hand, naturally following our direct definition of the space of closure (14. 4p . 
which assumes the a priori knowledge of the expansion coefficients {c^} for any v G L^, 
such extended linear functionals for every (3 act as follows: 

(5.15) {v,ip*p)^ = {v,ip*p) denote the expansion coefficient of v in (14.51) . 

In view of the performed construction of via closure of finite sums, it is not difficult to 
see that such continuous extended linear functionals are defined uniquely (by the density 
of finite sums in (14.41) ). and in L? are given by the standard integrals as in (13. lip and 
(I3.22p . Defining a countable and complete set $* of such extended linear functionals 
suffices for our applications. 



^"If X is a linear normed space, i is a linear manifold, and / is a linear continuous functional defined 
on L, then / can be extended to F on X and = 
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Thus, overall, in the sense of (15. 151) . the standard bi-orthonormality of the bases {ipp} 
and becomes trivial: 

(5.16) {iIji3,iIj*)* = {'ipi3,i'^) = ^is-y for any (multi-indices) /3 and 7, 

where (■, ■) is the usual duality product in L^(M) and 5/?^ is the Kronecker delta. Then, we 
treat similarly the adjoint linear functionals {w,iljp)^: for any w G 1/^,, where $ = {ipjs} is 
now a countable (complete) set of linear functionals, which have been properly extended 
from the adjoint space L'p, {M.) = /^^.(ffi), with p* = = p. 

According to the above analysis, we now show some convenient formal calculus with 
non-convergent integrals, which illustrate standard bi-orthonormality and other proper- 
ties. For instance, according to these formal integration techniques, we have: 

Proposition 5.3. 

= 1 for all /3>0. 

Proof. We start by looking at the scalar product, in the indefinite metric, defined by 

= - 1 i^fiiyWpi-y) dy- 

By our definitions of i'jsiy) (15.131) and ip^{y) (15. 5p . we substitute to find that 

r L2fc + lJ -| 

L j = l -I 

We now apply to ip^ the identity operator (D^)^^D^ = /, with a standard definition and 
construction of the inverse integral operator {D^)~^ = D~^, such that 

(^^,^;K = (^^,(D;^d,OV';)*- 

Integrating by parts, (3 times, we find that 

{^P,rp). = D^^F{y){D-^ D^^) ^ {{-y)f^ + o{y^)) dy 

= jr /(-I) [DS-'Fiy) D-/^^ [{-yf + o(/))] dy 
= - = ir /(-l)^ [nv) D'y {{-yf + o{y^))] dy 
= lj{-lfF{y){-lYf3\dy. 

One can see that, according to these formal calculus of integration by parts, we each time 
improve the convergence properties of the integrals involved, meaning using a distribu- 
tional treatment of those integrals as values of certain linear functionals as generalized 
functions (distributions). As customary, this corresponds to a regularization of divergent 
integrals. 

Hence, eventually, it follows that 

{ipp,rp)* = ! F{y)dy = l, 

so finally we arrive at a convergent integral, but not absolutely, since F is not Lebesgue 
measurable in M, as we have seen. □ 
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Proposition 5.4. In terms of the above formal calculus, 



(5.17) {i;p,r,)* = ^ for all /? ^ 7. 

Proof. The first part of tlie proof follows that of Proposition 15.31 First, consider the case 
where /3 > 7. After integration by parts, it can be seen that {ipp.ip*!)^ may be written as 

{i^.^r,). = (-i)"/^Kr^^?^(2/)^^?((-2/r + •••) ^y- 

However, since P > 'j, then it is known that D^^^—y)'^ = 0. So it follows that (15.170 holds. 

Now consider the case when f3 < 'j. In this case it can easily be seen that the above 
argument will not work. Rather than attempting to use a similar argument, we instead 
use another proof which encompasses both cases of /5 > 7 and /5 < 7. 

By the definitions of the linear operator B and the adjoint operator B*, we know that 



(5.1^ 



which defines the eigenvalue problems for these two operators. Taking the inner product 
of (15.181) . in the indefinite metric, with ip* and ipp, respectively, yields 

However from the definition of the adjoint operator B*, we know that 

Hence, if /? ^ 7 (i.e., ^ A^), (ISlTTj) follows. □ 

6. Semilinear dispersion PDEs: VSSs, numerics, and preliminaries 

6.1. Semilinear model. We now consider the odd-order problem, but now with a non- 
linear absorption. We look at the Cauchy problem for the semilinear odd-order equation 
(11.201) . for k = 1, 2, 3, . . ., with sufficiently good initial data m(x, 0) = uq{x). As we have 
already commented on, this non-integrable model is somehow connected to the KdV equa- 
tion (II. 4p . with the difference being that the extra operator — which corresponds 
to absorption, is simpler and of zero differential order. This kind of nonlinearity also 
allows us to avoid entering remarkable classes of integrable PDEs, which obey various 
extremely strong and quite specific properties that are illusive for more general nonlinear 
odd-order PDEs. 

Let us mention again, as in the linear case, even-order semilinear problems have been 
studied over recent years and are fairly well understood, unlike similar odd-order ones. 
The related generalized semilinear even-order model is given by 

(6.1) = -(-A)™?i± |u|P-V 

which includes both the cases of absorption and reaction. The lower order absorption 
case, with m = 1, corresponds to the heat equation with absorption given in (11.131) . 
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6.2. Similarity solutions of semilinear equations. As is customary in PDE theory, 
there exists a critical exponent for (11.201) given by 

(6.2) p = p, = i + ^\^^^ = 2k + 2. 

It can be called the critical Fujita exponent; see further comments below. As for the 
parabolic equation ( 16. ip . where 

Po = 1 + ^ 

(as in (16.21) . 2m stands for the order of the differential operator involved), the critical 
Fujita exponent characterizes parameter ranges of blow-up and non-blow-up solutions and 
changing of the stability of the trivial zero solutions for the PDEs under consideration. 
As usual, we consider self-similar solutions of the very singular type of (ll.20p 

u^{x,t) =r^f{y), y = x/t^, 

where / solves the ODE 

(6.3) (_i)^'+i/(2m) + _i_j/^ + _i_j_|j|P-ij = in M. 

Unlike the linear case, we cannot integrate to reduce the order of the ODE and hence 
we remain with an equation of order 2fc -|- 1. We see that (16.31) is a difficult higher-order 
equation and so begin with numerical results. 

6.3. Numerical results for the semilinear equation: k = 1. We look at similarity 
profiles of the semilinear equation (11.61) . for k = 1, 

(6.4) /" + i/'^ + _i_/_|/|P-i/ = in M. 

Once again we use the MatLab bvp4c solver to plot the profiles. This however gives us 
the same problem, as in Section 12.31 of solving an initial value problem using the BVP 
solver. We take / and /' to be zero, as boundary conditions at the left-hand boundary 
point of a fixed interval and / to be zero as the condition at the right-hand point. This 
corresponds well to a proper BVP setting for odd-order operators; see p^y as a guide. We 
then try to shift the point at the right boundary in order to coincide with a point where 
the oscillations go through zero and hence find the best profile. We look in particular for 
convergence and refiectional (about the axis / = 0) symmetry of the tail. 

Figures H] and O show two sets of profiles for p = 1.9 and p = 3.3. Whilst we cannot 
guarantee the accuracy of the size of the tail, the profile close to the origin, where max |/| 
is obtained, is very stable. We see that as p decreases, max |/| increases and this seems 
to justify the term "Very Singular" , with the mass concentrated close to y = . 

In the linear case set out in Section 12. 3[ due to the instability of the solutions, we had 
to use a "matching technique" in order to find reliable profiles. However, in the semilinear 
case, this is not necessary, as the profiles are stable, since we do not have a scaling group 
of solutions {c/}, due to the nonlinear term \f\^~^f- In fact, we find that whilst the 
tail of the profile may differ, given different boundary points, the rest of the structure is 
extremely stable and rarely changes. For, approximately. 



2 <p < 3, 
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we can find profiles easily, for almost every boundary point value we use and they only 
really differ from the tail. But for other values it is more difficult to find profiles, especially 
reliable ones. 

In particular, we look to see the behaviour of the solutions as p — 1^ and as p — 4~, 
which is the critical Fujita exponent (16. 2p for k = 1. 

Since we are using the boundary value solver of MatLab, we face numerous difficulties 
in obtaining the correct numerical results. As we are solving the initial value problem 
of the semilinear equation as a boundary value problem, difficulties arise in finding the 
correct right-hand end-point of the interval. 

As we do not know exactly at which points / and /' are both zero, we attempt to 
approximate the point by looking at various profiles and finding the best. If the far 
right boundary point is incorrect, then we produce artificial oscillations. We know which 
profiles are most likely to be false from the analysis we have done on odd-order linear 
PDEs. Since the semilinear equation can be naturally treated as the linear one with 
a perturbation — we expect similar behaviour. With a wrong boundary point, 
we can end up with non-symmetric tails, since we are forcing the oscillations through a 
specific point. Hence, the most important condition we look for in reliable profiles are 
that the oscillations are symmetric. 

We also look to ensure that the oscillations become symmetric as quickly as possible. 
This also follows from having forced oscillations due to wrong boundary points. Finally, 
we look for oscillations that decay as y increases. 

The solver that we use, bvpinit, requires as the name suggests an initial guess for the 
function. Whilst this guess does not necessarily have to be very accurate, occasionally a 
wrong guess can lead to wrong profiles being found. We have yet to find any evidence so 
far, to suggest that the initial guess is bad enough to affect the profiles though. However 
for less stable profiles where p ~ 2, the solutions blow-up and the initial guess for max |/| 
may not be accurate. The values needed for the initial guess for the maximum value 
are much larger than the maximum value found in the profile. Whilst the profiles look 
reliable, it is unclear as to whether the inaccurate initial guess affects the output. 

Thus, the conditions we look at to ensure the best possible profiles can be briefly 
summarised as follows: 

• Symmetry (reflectional, / t-^ — /) of tail for y ^ 1. 

• Symmetry of tail occurs as close to as possible. 

• Minimisation of symmetric tail. 

Finally, in Figure [6] we present the VSS profile for p = 3.99, which is sufficiently close 
to the bifurcation critical Fujita exponent Pq = 4; see Section [73] for derivation. 

6.4. Numerical results: k > 2. For k = 2, scaling (12.11) in the PDE yields the following 
equation for /: 

(6.5) ut = -u^.,^^-\u\P-^u =^ _/(5) + i/'^ + _i_/_|/|p-i/ = o in R. 
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(a) General view of f{y)- (b) "Tail" of the solution, for y ^ I 

Figure 4. A VSS profile of <^ for p = 1.9. 




-0.1 I ' ' ' ' ' ' 1 -0.04' ' ' ' ' 1 

50 100 150 200 250 300 275 280 285 290 295 300 

y y 

(a) General view of f{y) (b) "Tail" of the solution, ioi y ^ 1 

Figure 5. A VSS profile of for p = 3.3. 



In Figure [TJ we show the similarity solution of the ODE (16.51) for p = 4:, where the tail is 
not that oscillatory as used to be for /c = 1. Figure [8] shows VSS profiles for p = 3 and 
p = 5.9, which is slightly lower than the critical exponent Po = Q for k = 2. 

Finally, we present a few VSS profiles for the seventh-order equation for A; = 3: 

(6.6) Ut = u^,,,.,.-\u\P-'u ^ 1/'^ + _!_/_ |/|P-i/ = o in M. 

These are given in Figure [9] for p = 4, 5, and 6. Note that the first critical exponent is 
Po = 8 for k = 3, so these values are well below the po-bifurcation point; see Section [731 
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Self-similar profile for k=1 , p=3.99<Pjj=4 
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Figure 6. A VSS profile /(y) of ([63]) for p = 3.99 < po = 4. 



Self-similar profile for k=2, p=4 
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Figure 7. A VSS profile f{y) of ([63]) for p = 4 < po = 6. 
7. Linearized stability analysis and centre subspace behaviour: 

APPLICATIONS OF HERMITIAN SPECTRAL THEORY 

7.1. Linearized stability and critical (Fujita) exponent. We begin with some ap- 
plications of the spectral analysis obtained before, in the theory for LDEs. 

As in the linear case, we introduce the following similarity scaling in the semilinear 
equation (11.201) . with 

t) = (1 +t)"2fcTT-y(y^r), y = x/{l + t)'2ih^, r = ln(l + t). 
Our rescaled equation is then given by 



(7.1) 



32 



^— V — \v\^ 
p-1 I I 



Biv 
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Self-similar profile for k=1 , p=3<p =4 
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Figure 8. VSS profiles f{y) of ([63]) for p = 3 (a) and p = 5.9 < = 6 (b). 
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Figure 9. VSS profiles of the ODE 1^ for p = 4, 5, 6. 
Here, the linear operator Bi is defined by 

B, = B + dJ, where = ^ - ^ = 

and B is our canonical rescaled one (12. 2p . 

Lemma 7.1. For p > po = 2k + 2, zero is exponentially linearly stable for (17. ip in 



i7f+i(M)nL2(R). 

Proof. We look at the linearized problem of ( IT.ip about zero, 

f,- = Biw. 
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Then the spectrum is given 

by (t(Bi) = {di - for /> O}. Therefore, for / = 0, 

Xq = di = — < 0, for all p > 2k + 2. So, in our space, 

\\v{r)\\2k+i,p ^ e^°^ ^ as r ^ oo. 
Hence zero is linearly stable for p > pq. □ 

Lemma 7.2. For 1 < p < po, zero is exponentially linearly unstable for (17.1 1) , for small 
data vo e i/f+i nL2(M). 

Proof. For I = 0, Xq = di = — -^j^ > 0, so that 

||^(r)||2fc+i,p ~e^«" ^ +00. 

Hence zero is linearly unstable for p < pq. □ 

Thus, po = 2k + 2 plays almost a full role of the critical Fujita exponent for (11.201) . that 
mimics the standard one po = 1 + ^ (= 2m + 1 for = 1) for the semilinear parabolic 
equations (16.11) : see [13] for extra details. 

7.2. Centre subspace behaviour. Let us show that, at some critical values of p, there 
exists some quite special asymptotics of solutions. 

First critical exponent pq. We look again at the rescaled equation given by (17. ip . 
Let us look at the case where / = 0, with critical exponent p = Pq = 2k + 2. We check the 
behaviour close to the centre subspace of B, i.e., 

v{t) = co{t)iIjo + v^{t), 

where Vq is asymptotically small in comparison with the first term and is orthogonal to 
ipQ, i.e., {vq,iPI)^ = 0, for any k >1. Then, since Bipo = 0, we multiply the equation by 
t/'q = 1 to get the following leading term: 

co = -icor'+'co(i^or'+Vo+...,^^s)*- 

We need to assume that (see a justification below) 

7o = (|V'o|''+Vo,V'o)*>0. 

Note that analytically, proving that 70 7^ is difficult and even checking this numerically 
is also questionable. So, assuming that 70 > 0, 

\co\-^''^'^'f^=-lo + ... =^ -^|cor(2'=+^) = -7or + ... as r ^ 00. 
Finally, this yields the following rate of decay: 

Co(r) ^ [{2k + 1)70 r] for r > 1. 

Thus, the centre subspace behaviour of u{x,t) is given, as t — 00, by 

u{x,t) ^ (l + t)"5fcTT[(2fc + l)7oln(l + t)]"^^o(a;/(l + i)^), 

i.e., contains a typical extra logarithmic factor. A full justification of such a behaviour 
remains open (even existence of an invariant centre manifold is obscure). 

34 



Other critical exponents pi. Let us now look at the general case for /, with critical 
point p = pi = 1 + ^pf-, where we check the behaviour close to the one- dimensional kernel 
of B — A//, by setting 

V = ci{t)iPi + v^. 

Here, is small and orthogonal to ipi, as before. Then for Bipi = Xiipi, we take the inner 
product with ipi to get, for r 1, 

2fc+l 2fc+l 

c'l ^ -|q| '+1 ci{\^i\ '+1 ^pi,ipt)^. 
For convenience, let us assume, as usual, that 

, 2k+l 

Hence, assuming also that Q(r) > 0, we have, for r ^ 1, 

_2fc+2+i 2fc + l 

Hence, the stable subspace behaviour of u{x,t) as t — > oo, for all critical exponents 
p = Pi, is given by 

u{x, t) ^ (1 + t)"^ 7« ln(l + t)] (x/(l + t)^) . 

So, there exists a countable set of stable subspace behaviours governed by eigenf unctions, 
corresponding to the point spectrum cr(B) = { — 2j^^ ^ ^ 0}. 

7.3. Why 7/ > 0: numerics. For the above analysis to hold, it is assumed that 

, 2fc+l+i 

(7.2) ^l = m^T^^|Jl,^p;),>0, 

which is not an easy inequality to prove. Note that, for all / > 0, the eigenf unctions ipi{y) 
and ipi{y) are of the same parity, so that the indefinite metric (■, ■)=„ plays no role here 
and can be replaced by the standard L^-one (for odd Ts, with the minus sign). 
Let us first look at the case Z = 0. We need to show that 

70= (|V^o|''+Vo,V'o)* >0. 

However, in this case we have that ip^ = 1, so m. essence it is enough to prove that 

70 = /|^oP'+Vo>0. 

Attempting to prove this rigorously, is very difficult as well. However, numerically it can 
be shown, for A; = 1 at least, that this inequality is true. 

Certainly, we know that the first eigenfunction ipQ = F{y) can be found, for the lower- 
order case k = 1. Using the MatLab function trapz, which uses a trapezoidal method of 
integration, the integral can be approximated to 

(7.3) j |F|3F = 0.0300... > 0. 

It is noted that, since this numerical solution F, is not the fundamental kernel satisfying 
JF = 1, we scale the calculations such that this is true. So (17.31) holds for the true 
rescaled fundamental kernel. 
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Solving the problem (17. 2p for A; > 2 is more difficult as the shooting problem to find 
further eigenfunctions is not easy, but not impossible. Similarly, it is possible to construct 
solutions for all iph where / > 0. Again the problem of shooting is difficult, but not 
impossible to justify the inequality (17.21) . 

7.4. Bifurcation points and p-diagrams. Once again, we look at our semilinear equa- 
tion given in (16.31) . which we can write now as 

B/ + (ft - ^Fa)/ - l/rV = O^BJ- l/rV = 0, 
where Bi = B + dil is defined as before. Critical exponents {pi} occur when 

1 1 _ \ ; ^ „rU^^ 1 1 _ « 



di = — r - WTT = -^h i-e., when 



Pl-l 2k+l " ' 2fc+l 2A:+1' 

Therefore our critical exponents, pi, are given by 

P; = l + ^, / = 0,1,2,.... 

We can see from this that — > 1+ as / — > +oo. 

We look at p near these critical values, so that p ^ pi. We set e = pi — p, and then 

(7.4) (B~XiI)f + eaof = \fr'f + 0{e'), where ao = j^={^,)\ 

By classic bifurcation theory [56], our solution / can be given by 

f = CiJi + w^, 

where is orthogonal to ipi, 

Thus, taking the inner product of (17. 4p with ipi, we have that 
Since {ipi.ipl)^ = 1, we find 

(7-5) \Cr'= ^i^,i.-%,,^^). =i,{m'^^ ^here = 

For p = 1, we have that ni = {ipi,ipi)^ = 1, and so by continuity with respect to p, we 
must have that 

Ki > for all p ^ 1'^. 

Therefore, (17.51) indicates a countable number of subcritical pitchfork bifurcations. 

Figure [TO] shows a numerical calculation for the first branch of the bifurcation diagram, 
where we take / = and k = 1. Hence, in this case, the critical point is po = 4. 
During each iteration of the numerical program, the calculation uses the previous results 
to calculate the next step, thus improving the accuracy. The step size used here is 0.001, 
in the range of p = 1.7 to p = 3.3. Extending the range of values of p closer to the critical 
Po = 4: proves to be difficult. 

Figure [11] shows how the p-bifurcation diagram for the semilinear equation is expected 
to look like, if the branch is extended, given the numerics in Figure [TOl 

However the numerical analysis done in Section 16.31 has failed to provide firm evidence 
that Figure [11] shows how the p-branches behave, since we were not able to get reliable 
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Figure 10. Bifurcation branch for / = and k = 1, for the semihnear ODE (|6.3p . 
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Figure 11. A countable family of expected p-bifurcation branches for the ODE (j6.3p . 

numerics close to po = 4. As before, this is related to the extremal oscillatory behaviour 
of similarity profiles at the right-hand side, which does not allow us using standard nu- 
merical codes of continuation with respect to the parameter p, i.e., numerically construct 
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the so-called p-branches of solutions. In Figure [TTl we use the analytical evidence of 
such bifurcations from zero a.t p = pi, which also requires extra difficult mathematical 
justification. 
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Appendix A. "Radiation conditions" in spectral {B, B*}-theory 

Here, we clarify the "radiation-type conditions" posed at infinity in the domains of operators 
B and B*, which allow the operator pair {B, B*} to have purely discrete spectra { — 2k+i i ^ ^ 0} 
already detected by eigenfunction expansion of the corresponding semigroups. This serves as a 
description of new properties of the spaces of closure -C/p(M) and L^*(M). 

A.l. Domain T){&) of the linear operator B. To this end, in order for our linear ODE 
to be "well-posed" (i.e., with a proper number of boundary conditions at infinity), both in a 
mathematical sense, as well as a physical sense, we look for conditions, which must be satisfied 
for the eigenvalue equation: 

BVz(y) = AjVi(y) in 

This can be rewritten as 

(A.l) {-if+^,i^f'^'\y) + ^,Uy) + 2FH yV'Ky) = ^My)- 

Since the order of this ODE is 2k -\- 1, it is natural that there must also be 2A; -|- 1 boundary 
conditions placed, as classic theory of ordinary differential operators suggests; see Naimark's 
monograph [46j . 
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Eliminating exponentially growing bundles . First consider the problem, as y ^ +00. 
Attempting to balance leading order terms in (jA.ip . leads to 

(A.2) + ^^y^l,[{y) ~ 0. 

As y — > +00, we have that 

2fc+l 

(A.3) V'/(y)~e^^^, 6gC, 6^0, 

hence substituting this into the above equation, yields 

It can be seen that we end up with two cases, dependent on the value of the parameter fc, 
which will determine the sign of 6^'^ and therefore its roots. Hence, for now, we ignore the term 
( 2fc+i ) 2k+i ^'^'-^ j^^* \ook at the value of (—1)'^, assuming that b = ( 2^+1 ) 2k+i 

When k is even: for even values of k, it is noted that b'^^ = 1. Hence there must 2k roots for 

11 1 7 27rmi 

0, which are given by Om = e 2fc for m = 0,1, ... ,2k — 1. 

When k is odd: for odd values of k, there is now a negative sign, such that b'^^ = —1. Similarly, 

^ (7r + 27rm)i 

as before, we derive 2k roots, where now b^n — e for ui — 0, 1, ... , 2k — 1. 

For the problem to be well posed on the space L^, it is important to look for roots such that 
there is exponential decay, rather than growth. In other words, it must satisfy the condition 
p = e~"^", as y ^ +00. Hence we need to eliminate any roots such that Re^^ > 0. It is first 
noted that equality, Re 6m = 0, occurs when 

for even k, 

for odd k, 

with the same applying for For Re bm = 0, we must have that 

{I for even k, f ^ for even k, 

^ for odd A; for odd /fe. 

Hence, in order to eliminate roots which give rise to exponential growth, the conditions placed 
must be such that we do not include roots such that 




{m < -I and m > ^ for even k, 
m < and m > for odd k. 

Therefore, by taking the weight 

2fc+l 

(A. 4) p{y) = Q~^y^^ ^ with any sufficiently small a > 

we eliminate all exponentially growing oscillatory bundles. A sharp bound on admissible a > 
will be derived. Thus, we have k — 1 conditions placed as y ^ +00. 

Similarly, we can do the same analysis for y— >— cxd. Asy— >— cxdIu (IA.2p . we have that 

2fc+l 

(A.5) ^(y)~e^(-^)- ^ t^'^ ~ ("1^ ^ (ife) 1^- 
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As expected, this only differs from the y +00 case by the opposite sign. This leads to 

{(7r + 27rm)i 
e 2fe , for even fc, 
jrmi 
e k ^ for odd k, for m = 0, 1, . . . , 2A; — 1. 

As before, we do not want roots such that Re bm > 0, where in this case 

m < and m > '^^2^ f*^^ even k, 

m < I and m > ^ for odd k. 

Hence, for the weight (j3.17p . we eliminate these conditions. This leads to a further k condi- 
tions, which are placed at y ^ —00. All these correspond to eliminating exponentially growing 
asymptotic bundles that is not enough as we explain below. 

Domain P(B): radiation condition. As stated before, we look for 2k + I conditions to 
be posed onto the problem. Hence there are two more conditions needed. These conditions 
at y — > ±00 are known as radiation conditions. In classic problems of quantum mechanics, 
acoustics, and physics, the general idea behind radiation conditions is that energy sources must 
exactly be that and not sinks of energy. Hence all energy must be radiated from a point 
and scatter to infinity. We refer to the book by Sommerfeld [50], who (in 1912) first proposed 
radiation conditions for the Helmholtz equation, and to [5^ for recent applications and references. 

In our problem, the radiation conditions are rather tricky and have almost nothing to do with 
the classic ones. We recall that we identify those just for convenience (to verify the domain 
of B and B*), since the eigenfunction expansions of the semigroups, as the main tool of our 
asymptotic analysis, automatically include the necessary two conditions at infinity. 

Remark. As is shown in |21j, for the linear rescaled Schrodinger operator in L^(]R^), p = 
e~"l^l", a = (a standard setting), 

B = -i (-A)- + ^ y . V + ^ I, B* = -i (-A)- - 2k y • V, 

the absence of a "radiation condition" , accepting "oscillatory functions" at infinity only, will con- 
sequence in the fact that, in addition to the standard required spectrum o"(B) = {— > 0}, 
there appears another symmetric real one {+2^, \P\ > 0} with generalized Hermite polynomial 
eigenfunctions constructed as in (j5.5p . Fortunately, this is not the case for our operators (j2.2p . 
since the polynomials for the weight in (j3.17p . 

The origin of the radiation condition for B is as follows: Let us now balance all lower-order 
terms in the eigenvalue problem ()A.ip . so that 

By integration we can easily see that 

(A.6) Vi(y)~^y^''+'^''-\ 

for some constant A. Of course, this corresponds to the obvious root 6 = in the exponential 
expansions ()A.3P and (jA.Sp . 
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We note that (jA.6p . is a "rational" function, unlike the exponentially oscillatory bundles in 
(1A.3P and (jA.Sh . For y —oo, we know that, for any A G C, rational solutions such as (IA.6P 
do not belong to the space with the exponentially growing weight (j3.17p . 

Thus, overall, we conclude as follows: 

(A. 7) at y = — oo, the proper weight (j3.17p generates k + 1 conditions. 

So, this is a usual and a standard situation, so that the singular point y = — oo does not require 
any radiation- type condition. This is not the case for the "oscillatory" end-point y = +00. 

For y — > -|-oo, consider all complex "eigenvalues". A; G C, such that A; = P + iQ, for some 
PjQgM. From (lA.6p . ^i{y) may now be given by 

(A.8) ^z(2/)~/+"^, 
for P, Q G M. Hence we see that 

(A.9) i^iiy) ~ yV^''"S' ~ /(cos(61ny) + isin(61ny)), 

as y — > -|-oo. However, we know from the asymptotic analysis that the behaviour of proper 
eigenfunctions is different and given by a distinct type of higher oscillatory functions: 

/ r \ / 2k+i 

My) cos(y 2k ), 

which obviously gives a stronger oscillatory behaviour than a pure cos (6 In y) in (lA.Op . However, 
(jA.OP admits weaker oscillatory behaviour and so we must place a condition to eliminate this 
behaviour. We recall the all the proper eigenfunctions being given by the generating formula 
''Ppiv) = '^0fD^F{y), do not contain the bundle ()A.6p . since the fundamental rescaled kernel F{y) 
does not by the known divergence of the operator B (the equation for F has been integrated 
once with the zero constant of integration that eliminated any trace of (IA.6P ): see below. 

Thus, the generalized radiation condition, that is necessary for the proper discreteness of the 
spectrum of B in L'i, can formally be formulated as follows: 



(A. 10) for the eigenvalue equation (lA.ip . the bundle (jA.Sp at y = +00 is absent. 



Then, as usual, the domain By the construction of the space of closure in Section 14. 3[ the 
condition (|A.10p is assumed to be included, though we do not prove this here; cf [21] for linear 
rescaled Schrodinger operators. 

Actually, it is easy to see that all our eigenfunctions {ipi{y)} satisfy the above radiation 
condition. Indeed, F{y) = V-'o(y) does satisfy this, by integrating once, where we have that 

for some constant of integration C. The last term precisely shows that such a rational behaviour 
is absent, since we have that 

F{y) ~ 

which implies that C = 0. Then each eigenfunction 'ipi{y) = F^^'^y) for all I > 0, also satisfies 
this condition. 
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A. 2. Radiation condition in the domain D(B*) of adjoint operator B*. We apply the 
same analysis to the adjoint operator B*, where the eigenvalue problem is now 

(A.ll) i-l)''+^^Pt^^'+'\y) -^^yr'iy) = Xiriy)- 

Similarly, as before, we look at exponential bundles with b ^ 0: 

2fc+l 2fc+l 

ip*{y) ~ e^^ as y — > +00, and V'*(y) ~ 6^'^' as y — > —00. 
Balancing terms 

(A.12) = iFTT y^'iy) + 

it can easily be seen that the case for the adjoint operator B*, only differs from the linear 
operator B, with respect to a change of sign. Hence it is seen that, as y — > +00 and y — > —00 
respectively, for m = 0, 1, . . . , 2A; — 1, 

{ (7r + 27rm)i / nmi 

e 2k , for even k, ^ _ I e fc , for even k, 

vrmi and bm = S (,r + 27rm)i 
e , tor odd k, , tor odd k. 

It is noted that we look at the problem in the space L^, , with weight 

2fc+l 

p*{y) = e~"l^l , a > is small enough. 
Hence, for this weight, for y +00, we eliminate the roots such that 

j m < and m > for even k, 



and for y —00 



m < ^ and m > for odd A;, 



m < I and m > ^ for even /c, 
m < and m > for odd A;. 



These give 2k — I conditions in total. 

Whilst before balancing the rest of the terms will lead to our radiation conditions, this is not 
possible to do this in similar lines (as for B) in the case for B*. Note that, on integration, there 
holds: 

for some constant Aq. However, this behaviour is perfectly acceptable in L^, , and it is satisfied by 
the polynomial adjoint eigenfunctions {V'f (?/)}) which are our generalized Hermite polynomials. 

Instead we look for the remaining two conditions for the problem as y — > —00. In light of the 
definition of B* , we assume to have another type of radiation condition, which for convenience 
we denote as the "Adjoint Radiation Condition" , such that we 



exclude two bundles with "maximal" oscillatory components at y = —00. 



Again, this extra radiation condition makes the total number of conditions to be equal to 2A; + 1, 
which is the differential order of the operator B*, so that the eigenvalue problem becomes 
algebraically consistent [l6]. In other words, we have an algebraic inhomogeneous system of 
2A; + 1 equations with analytic coefficients with 2k + 1 unknowns. Such systems do not have 
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more than a countable set of solutions, which are eigenvalues of B*, which is defined in such a 

way. 

Hence, we now restrict those roots such that Rc a„i = 0. Hence the conditions, as y — ^ +00 
are now given by the following distribution of the acceptable coefficients {5^}- 

{m < I and m > ^ for even k, 
m < ^ and m > for odd k. 

A.3. Calculations for the weights, p{y) and p*{y). We now determine the sharp distance 

between the principle root such that Rc b„i^ = and the previous root bmc-ii where 
Re 

bruc-i > 0- III doing so we may find our weight p{y) such that it cuts off all unwanted 
roots, for which Re 6 > 0. 

First consider the case as y ^ +00. For Re = 0, it is known that there is a root here and 
this is given by 

, , . I I for even k, 

bmc = dkh where dk = 2K;(2FFTj "^c = 



for odd k. 



Hence, the root rric — 1 is given by 

rric-l 

This yields 



for even k, 
for odd k. 



bmc-i = {(^os^^-^ + ism^^-^)dk, 
for all k. Hence the distance between the two roots (in the real axis) is 

d = dk cos -^^^ > 0. 

For y —00, we do not have any roots Re bm^ = 0, for any A;. However we look at the 
distance between the real axis and the next root such that Re 6 > 0. 

In this case we find the distance between the real axis and roots such that 

for all k. Hence we have that 

'4[cos(L^Jf + t) + isin(L^Jf + i)] for even k, 
.4[cos(Yf)+isin(^f)] for odd A;. 

This gives the distance between this root and the real axis 

^ _ f 4 cos( [^J f + 5^) for even A;, 
[dfc cos(^|) for odd A;. 

This characterises our weighted space LVM.), with the exponential weight 



■'mc—l 



p{y) 



2k + l 

^a\y\^ for y<-l 

2k+l 

^-ay 2k for y > 1, 
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where a G (0, 2d). 

Similarly we can find the weight p*{y). We note that the calculations are exactly the same as 
in the case for p{y)^ except a difference in sign when calculating the roots. This leads to 



d 



as y — > +00 and 



4 cos([ V-J f + ^) for even k, 
dk cos(^^l) for odd k. 



dk cos (^2^f ) for even fc, 

k-2 



[4 cos([^^Jf + ^) for odd A:, 
as 2/ — > — 00. Here the weight p*{y) is defined as in (|5.4p . 



Appendix B. Estimates on the fundamental kernel and majorizing 

operator 

Here, we develop some ideas concerning majorizing-comparison issues for the odd-order opera- 
tors involved, which literally do not have any connections with the standard Maximum Principle. 

B.l. Estimates of the rescaled kernel. Recalling our estimate given by (|2.1ip . we now look 
to estimate our rescaled fundamental kernel F{y) by 

(B.l) \F{y)\ <DF{y), where > and jF{y) = l. 

Here D \s a. normalisation constant, obviously satisfying D > 1. There exists infinitely many 
functions which satisfy (jB.ip . but since our kernel F is changing sign, it is not possible to find 
an optimal analytic function 

Fo^t{y)=u:i\F{y% 

where wi > is a normalisation constant, such that 

/F„p, = l =^ u, = {J\F\)-^>l. 

However, we can find an analytical approximation of -Fopt such that (jB.ip is satisfied, but is 
non-optimal. One such function is given by 

A sketch of the function F^, is shown by Figure IB. II and comparison with the numerics in Section 
12.31 shows how the function may be an upper bound for F{y). 

By using the asymptotic analysis set out in Section 12.21 we can also find that all higher order 
derivatives are estimated by 

I —U !U 

B.2) \T^^F{y)\ < 



^I3p^l^e-''\y\°' for y < -1, 



c^y2k for y > 1, 



where c is dependent on k only. 

46 



F 


\x F ~ y-(2fc-l)/4fc 





y 



Figure 12. Majorizing kernel 



B.2. Majorizing kernel and spectral properties. We now introduce the positive majorizing 
"fundamental solution" 

with the majorizing F relative to the kernel F of 6(3;, t) in (12. ip . We then end up with a formal 
integral evolution equation written in the standard form 

(B.3) ut = A{t)u, 

for some linear operator A(t), with the "fundamental solution" b{x,t). This formal non- 
autonomous (in time t) evolution equation is understood in the sense that the Cauchy problem 
for ()B.3P with initial data 

(B.4) u{x,0) = uoix)>0 in M, 

is given by the convolution 

(B.5) u{x, t) = M{t)uoix) = b{t) *uo = t~2ih J F{{x - z)t~'2m:)uo{z)dz. 

It can be seen from this, that for general kernels F, majorizing semigroups do not exist, so 

equation (|B.3p does not admit translation in time. This defines the corresponding majorizing 

integral equation. For higher-order parabolic (poly-harmonic) equation, the idea of majorizing 

integral operators was introduced and applied in blow-up studied in ^23j . 

As before, let us now introduce rescaled variables 

1 _ 1 
u{x,t)=t 2k+iyr,(^y^r), fj = xt 2fc+i ^ r = lni. 

Hence from convolution in (jB.SP 

w{y, t) = / F{y — ze~ ^fc+i )uo{z) dz. 
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Using Taylor's power series we formally obtain 

(B.6) F{y-ze-^)=Z e'^^ ^-^D^F{y)z^ ^ ^ e'^^ ^ Mv)^^^ 

where Vi/3(y) = ^^D^F{y). 

The convergence of ()B.6P on bounded intervals is guaranteed by the estimates of DyF{y) given 
in ()B.2p . A proper convergence in can be also included as above. 
The solution can then be represented by 

w{y,T) = ^2^'"^ Mi3{uo)'4>f3{y), A/3 = -2^, where M/3(uo) = -4sj f z^^ uo{z) dz 

(13) ^' M 

are the corresponding moments of the initial data. 
Proposition B.l. There exists some formal operator B, such that 

Wr = Bu; 

and this induces the majorizing semigroup {e^"^}, formulated by the rescaled variables. 

It follows from (j3.1Up that B has point spectrum given by o"p(B) = {A/3} and corresponding 
eigenfunctions are thus given by 4> = {tpp}- With a standard definition of the space of ^>-closure 
Lp (see Section 1331) , this becomes the only spectrum of B. 

In order to find a full expansion form of the majorizing semigroup, we once again perform 
another rescaling given by 

u= {l + ty'^^w, y = + t)"2fcTT, r = In (1 + t) : M+ ^ M+. 

Then rescaling the convolution we obtain 

w{y,T) = e^^uo = (1 - e"^)~2feTT J F{{y - ze~^){l - e~'^)~'2^)woiz) dz. 

Using Taylor expansions, we find that the solution is given by 

wiy,r) = {l-e~-)-^ Z ^-^'^ D-FiO)j^ y-^^ie- - 1)'^^ 

x(l - e~^y^^ J z^'wo{z) dz. 

R 

This again shows the discrete spectrum, as well as some traces of generalized Hermite polyno- 
mials {ip*p}. 

The adjoint operator B* with such polynomial eigenfunctions {ip'^} occurs if we use the blow- 
up scaling 

_ 1 

u{x,t) = w{y,T), y = x{\-t) ^k+i ^ ^ = 

We thus obtain 

?D(y, r) = (1 - e~^) / F((ye" sfft - z){1 - e~'^)~^i^)wo{z) dz. 

R 
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Using Taylor expansion yields 



w{y,T) = (1 - e ^) 




(-1)' 



I3\u\ 



13 




1 



y'^ ^ j z'^wo{z) dz. 



iu~i3y. 



As above, further expanding of exponential terms here leads to the eigenfunction expansion, 
which determines finite generalized Hermite polynomials. This representation is rather technical 
and we do not present and analyze it since do not aim using those polynomials in what follows. 
Our main application of the majorizing operators is as follows: 

B.3. Comparison with majorizing problem. We look to see how the majorizing kernel 
relates to our real solution u{x,t). By looking at the convolution ()3.ip for the linear PDE, we 
can easily see that 



We introduce the majorizing kernel, since the structural behaviour of self-similar solutions 
of the majorizing evolution equation describes essential features of the solutions to the original 
PDE. Since we have estimated our solution, it can be possible to find properties of the real 
solution, which otherwise would be more difficult to achieve. 

The linear semigroup for B is not order-preserving, since the kernel F is oscillatory. Therefore 
we must have D > 1, which gives the order deficiency of the linear operator B and of the 
linear convolution operator M{t), given in (|B.5p . The actual defect is, in fact, characterized by 
D — 1 > 0. The linear semigroup for B would only be order-preserving if F did not change sign, 
since then we would have that F = F and so we would have D = 1 and so the defect would be 
zero. As we have shown, this is not possible for any A; > 1. 
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Hx,t)\ < \b{t)\ * \uo{x)\. 
Now looking at our estimate of the majorizing kernel, b{x,t), we have that 

\u{x, t)\ < Db{x, t) * \uo{x)\ < b{t) * uo{x), 



where we have to assume the following inequality for initial data: 
(B.7) D\uo{x)\ < uo{x) in M. 



Proposition B.2. // ([Ri]) and (|R7l) hold then 

\u{x,t)\ < u{x,t) in Rx M+. 
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